Extra force in f(R) modified theories of gravity 
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The equation of motion for massive particles in f(R) modified theories of gravity is derived. By 
considering an explicit coupling between an arbitrary function of the scalar curvature, R, and the 
Lagrangian density of matter, it is shown that an extra force arises. This extra force is orthogonal 
to the four-velocity and the corresponding acceleration law is obtained in the weak field limit. 
Connections with MOND and with the Pioneer anomaly are further discussed. 
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I. INTRODUCTION 

Higher-order curvature theories of gravity have re- 
cently received a great deal of attention in connection 
with the possibility of giving rise to cosmological models 
in the context of which one can address the issue of the 
accelerated expansion of the universe without the need 
of ad hoc scalar fields [3, S, 1 0, S @. These theo- 
ries involve corrections to the Einstcin-Hilbert action by 
considering a nonlinear function of the curvature scalar, 
f(R). Earlier interest in f(R) theories was motivated by 
inflationary scenarios as for instance, in the Starobinsky 
model, where f(R) — R — A + aR 2 was considered [7(- 
Other motivations include the search for wormhole-type 
solutions Q. In these studies, different approaches are 
used throughout the literature. These include the metric 
formalism, where the action is varied with respect to the 
metric; the Palatini formalism, where the metric and the 
connections are treated as separate variables; and the 
metric-affine formalism, which generalizes the Palatini 
variation, where the matter part of the action depends 
and is varied with respect to the connection 0]. 

Recently, it has been argued that most models pro- 
posed so far in the metric formalism violate weak-field 
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solar system constraints (To[, although viable models do 
exist Furthermore, it has been argued that higher 

order gravity may ex plai n the flatness of the rotation 
curves of galaxies fl2l. Il3l. [T3|. For instance, in the con- 
text of f(R) = foR n theories, the obtained gravitational 
potential is shown to differ from the Newtonian one due 
to the appearance of a repulsive term that increases with 
the distance from the center. The rotation curves of our 
Galaxy were studied, and compared with the observed 
data, so to assess the viability of these theories and to 
estimate the typical length scale of the correction. It was 
shown at first approximation, where spherically symmet- 
ric and thin disk mass distributions were considered, that 
a good agreement with data can be obtained with just 
the stellar disk and the interstellar gas. 

In this paper we aim to derive the equation of mo- 
tion for massive particles in a class of generalized grav- 
itational models in which the Lagrangian of the gravi- 
tational field is an arbitrary function of the curvature 
scalar. The study of the equation of motion is of funda- 
mental importance for the understanding of the structure 
and properties of gravitational theories. One of the most 
effective ways to test gravitational theories is by match- 
ing their predictions with the motion of real objects. For 
this purpose, we point out that the covariant conserva- 
tion equation for a symmetric energy- momentum tensor, 
corresponding to matter is not, in general, conformally 
invariant (l5j . One is led to relax the covariant conser- 
vation of the matter energy-momentum by considering 
a coupling between the matter Lagrangian and an arbi- 
trary function of the curvature scalar. It is interesting to 
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note that nonlinear couplings of matter with gravity were 
analyzed in the context of the accelerated expansion of 
the Universe [IB], and in the study of the cosmological 
constant problem [l^|. This is reminiscent of the situa- 
tion in scalar-tensor theories of gravity and also arises in 
string theory. Non-minimal couplings have also been ex- 
tensively considered in the literature, namely, between a 
scalar field and matter, including baryons and dark mat- 
ter [3 M, [H M, Hi SI HI- These couplings imply 
the violation of the equivalence principle, which is hi ghly 
constrained by solar system experimental tests [24l. |25j. 
However, it has been recently reported, from data of the 
Abell Cluster A586, that interaction of dark matter and 
dark ener gy does imply the violation of the equivalence 
principle [26|]. Notice that the violation of the equiv- 
alence principle is also found as a low-energy feature of 
some compactified version of higher-dimensional theories. 

In what follows, by considering a coupling between 
a function of the curvature scalar and the matter La- 
grangian, we show that in higher-order curvature theories 
of gravity the equation of motion of massive particles is 
non-geodesic. Thus, as shall be shown in Section [Til the 
equation describing the trajectory of the particle has an 
extra-force term, which is orthogonal to its four-velocity. 
The non-geodesic nature of motion is a distinct feature of 
these f(R) theories, found also in the context of a scalar 
field model with a suitable potential and proposed (27| as 
a solution for the Pioneer anomaly problem [28| . We shall 
discuss this question in Section IlIII Furthermore, it will 
also be addressed how our approach can be regarded as a 
covariant realization of a Modified Newtonian Dynamics 
(MOND) [2^| , even though free from the arbitrariness of 
its phenomenological version and somewhat simpler than 
its more recent realization which involves besides grav- 
ity, a vector and two scalar fields, the so-called TeVeS 
approach [30j ] . 



II. THE EQUATION OF MOTION IN f(R) 
GRAVITATIONAL THEORIES 



The action for the modified theories of gravity consid- 
ered in this work takes the following form 



S 



f 1 (R) + [l + Xf 2 (R)]C r , 



-gd 4 x, (l) 



where fi(R) (with i = 1,2) are arbitrary functions of the 
Ricci scalar R and C m is the Lagrangian density corre- 
sponding to matter. Note that the strength of the in- 
teraction between f 2 (R) and the matter Lagrangian is 
characterized by a coupling constant A. Analogous non- 
linear gravitational couplings with a matter Lagrangian 
were also considered in the context of proposals to ad- 
dress the cosmic accelerated expansion [lg], and in the 
analysis of the cosmological constant problem [TtJ ■ 
Varying the action with respect to the metric g^ yields 



the field equations, given by 

F x {R)R pv - ^fi(R)g„v - W l ^uF 1 {R) + g liV UF 1 {R) 
= -2XF 2 {R)£ m R pv + 2A(V A1 V„ - g^U)C m F 2 {R) 

+ [l + \f 2 (R)]TW , (2) 

where we have denoted Fi(R) = fl(R), and the prime 
represents the derivative with respect to the scalar cur- 
vature. The matter energy-momentum tensor is defined 
as 



rp(m 



S(g^) 



(3) 



Now, taking into account the covariant derivative of 
the field Eqs. ©, the Bianchi identities, VG M „ = 0, 
and the identity 



(□V, - V v D)Fi = R^ WFi 



one finally deduces the relationship 
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l + A/2 



g^c r. 



rp(m 
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Thus, the coupling between the matter and the higher 
derivative curvature terms describes an exchange of en- 
ergy and momentum between both. Analogous couplings 
arise after a conformal transformation in the context of 
scalar-tensor theories of gravity, and also in string theory. 
In the absence of the coupling, one verifies the conser- 
vation of the energy- momentum tensor (3lj . which can 
also be verified from the diffeomor phi sm invariance of 
the matter part of the action [U, I32l.l33|] . Note that from 
Eq. (0, the conservation of the energy-momentum ten- 
sor is also verified if f 2 (R) is a constant or the matter 
Lagrangian is not an explicit function of the metric. 

In order to test the motion in our model, we consider 
for the energy- momentum tensor of matter a perfect fluid 



rp(rn) 



(e + p) u^u v - pg^ 



(6) 



where e is the overall energy density and p, the pressure, 
respectively. The four- velocity, satisfies the condi- 

0. We also introduce the 



tions u^u^ 



1 and u^u 



projection operator = g^\ — u^u\ from which one 
obtains h^xu^ = 0. 

By contracting Eq. ([5]) with the projection operator 
, one deduces the following expression 

(e + p)SiiAU v V v u"-(V v p)(aX-u"u A ) 
XF 2 
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(C m + p) <y v R) (5 V X - u v u x ) = . (7) 



Finally, contraction with g aX gives rise to the equation 
of motion for a fluid clement 



Du a du a 



ds ds 



(8) 



where we have introduced the space-time connection r" , 
which is expressed in terms of the Christoffel symbols 
constructed from the metric, and where 



e +p 



\F 9 



1 + \fr 



(C m + p) V V R + V„p 



h au . (9) 



As one can immediately verify, the extra force f a is 
orthogonal to the four-velocity of the particle, 



= 



(10) 



which can be seen directly from the properties of the 
projection operator. This is consistent with the usual 
interpretation of the force, according to which only the 
component of the four-force that is orthogonal to the 
particle's four-velocity can influence its trajectory. 

Notice that massless particles do follow geodesies and 
therefore for them f a = 0. 



where C is an arbitrary vector perpendicular to the vec- 
tor /. In the following, we assume for simplicity, that 
C = 0. 

The mathematical consistency of Eq. (fl~3|) requires 
that /'fl^O, that is, vectors / and a cannot be orthog- 
onal to each other. We consider that both vectors are 
parallel. Therefore, we can represent the gravitational 
acceleration of a particle in the presence of an extra force 
as 



-• 1 I 1 2 c2\ a 

ajv = r {a - a N - f ) — 
2 ja 



(14) 



In the limit of very small gravitational accelerations 
a^r <c a, we obtain the relation 



(15) 



If one denotes 



III. THE ACCELERATION LAW IN f(R) 
GRAVITY 

To derive the acceleration law in the f(R) gravity, we 
start by assuming that the motion of a particle in a space- 
time with metric g^ v is given by Eq. ([8]) . The presence of 
the extra force f a implies that the motion of the particle 
is non-geodesic. For f a = we recover the geodesic 
equation of motion. The usual gravitational effects, due 
to the presence of an arbitrary mass distribution, are 
assumed to be contained in the term = r™ u^u" . In 
three dimensions and in the Newtonian limit, Eq. ([8]) 
can be formally represented as a three- vector equation of 
the form 



a = fljv + / 



(11) 



where a is the total acceleration of the particle, a/v is 
the gravitational acceleration and / is the acceleration 
(per unit mass) due to the presence of the extra force. 
If / = 0, the equation of motion is the usual Newtonian 
one, a — Un-i which for a point-like mass distribution is 
given by a = —GMr/r 3 . 

Taking the square of Eq. (fTTj) one obtains 



/ • 4 



1 / 2 
2<° 



a 2 N 



(12) 



where the dot stands for the three-dimensional scalar 
product. Equation ([12]) can be interpreted as a gen- 
eral relation which expresses the unknown vector ajv as 
a function of the total acceleration a, the extra force / 



and the magnitudes a 



l N 



and f 2 . From Eq. flTJ) one 



can express the vector <2/v, as one can easily verify, in the 
form 



1 

CLE 



1 

2f 



then Eq. (|15p can be immediately written as 



a N 



a E 



-a , 



(16) 



(17) 



which has a striking resemblance with the equation put 
forward phenomenologically in the so-called MOND ap- 
proach [291 ]. It then follows that a w y/oEON, and since 
ax = GM/r 2 , then a w y 'cleGM fr — v 2 g /r, where 
Vtg is the rotation velocity 
influence of a central force. 



of the particle under the 
Therefore, it follows that 
vf g — > v 2 ^ = y/asGM, from which arises the Tully-Fishcr 
relation L ~ as — cleGM, where L is the lumi- 
nosity that is assumed to be proportional to the mass 

Notice however, that in the framework of f(R) grav- 
ity, is not a universal constant as it depends on local 
curvature features. This might explain why it is some- 
what difficult to match the whole galactic phenomenol- 
ogy within the framework of MOND (see e.g. [34[ for 
a critical assessment). Nevertheless, this feature of our 
model opens up quite interesting possibilities as we will 
see next. 

Indeed, in general, the definition of a E , Eq. (|16p . al- 
lows one to formally represent the extra force as a func- 
tion of a and a e , that is 



, , , (18) 

O-E \ a E / \ a E , 

Hence, through Eq. (fl"6|) . Eq. (fT4|) can be rewritten as 
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where 

KiHfeH^v^+IP <20) 

Upon substitution into Eq. (|16p . we verify a 2 — v^/r 2 = 
ciEGM/r 2 , which yields for a^: 

aE = ki +2f - (21) 

Suppose now that / ~ GMa/r, where a is a constant, 
then in the large r limit, when / — * 0, cle ~ ot 2 is a 
constant, whose numerical value is determined by the 
physical properties of the extra force. 

Given the environmental nature of this extra force, 
only phenomenology can guide us in its identification. 
In the galactic context, it seems natural to identify 
with the ao = 10~ 10 m/s 2 , the threshold acceleration 
of MOND. On the other hand, in the solar system and 
its neighborhood, the Pioneer anomaly, whether a real 
effect unrelated with systematic effects (see e.g. [35[ 
and references therein), suggests that as = apio = 
(8.5 ± 1.3) x 10~ 10 m/s 2 . Interestingly, our approach al- 
lows for a unified explanation for these two problems and 
can account for the fact that the characteristic accelera- 
tion of each class of observations is somewhat different. 



IV. DISCUSSION AND CONCLUSIONS 

In this work we have studied a class of generalized grav- 
itational models, in which the Lagrangian density of the 
gravitational sector is an arbitrary function of the scalar 
curvature and an explicit coupling between the scalar 
curvature term and the matter Lagrangian density. In- 
terestingly, we have found that the equation of motion of 



massive particles is non-geodesic. Therefore, the equa- 
tion describing the trajectory of particles exhibits a term 
representing an extra force, which is orthogonal to its 
four-velocity. We have also shown that our models have 
similar features with the phenomcnological approach of 
MOND, providing an alternative formulation to this pro- 
posal without the need of introducing, such as in the 
TeVeS proposal, extra fields besides the metric and the 
scalar curvature, which now plays the role of an addi- 
tional scalar field. Furthermore, we have shown that 
the extra force is consistent with the so-called Pioneer 
anomaly. A distinct feature of the formalism outlined in 
this work is that it allows to establish a connection be- 
tween the problem of the rotation curve of galaxies, via a 
solution somewhat similar to the one put forward in the 
context of MOND, and the Pioneer anomaly, even though 
the characteristic acceleration of these two classes of ob- 
servation is somewhat different, about 10 -10 m/s 2 . Cer- 
tainly, a more detailed study of the solar system implica- 
tions of our model via the parametrized post-Newtonian 
analysis remains still to be performed, but it will be con- 
sidered elsewhere. 
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